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ABSTRACT 



Context. SN2011fe was detected by the Palomar Transient Factory on August 24th 2011 in M101 a few hours after the 
explosion. From the early optical spectra it was immediately realized that it was a Type la supernova thus making this 
, event the brightest one discovered in the last twenty years. 

Aims. The distance of the event offered the rare opportunity to perform a detailed observation with the instruments 
on board of INTEGRAL to detect the 7-ray emission expected from the decay chains of 56 Ni. The observations were 
performed in two runs, one before and around the optical maximum, aimed to detect the early emission from the decay 
of 56 Ni and another after this maximum aimed to detect the emission of 56 Co. 

Methods. The observations performed with the instruments on board of INTEGRAL (SPI, IBIS/ISGRI, JEMX and 
OMC) have been analyzed and compared with the existing models of 7-ray emission from such kind of supernovae. In 
this paper, the analysis of the 7-ray emission has been restricted to the first epoch. 

Results. Both, SPI and IBIS/ISGRI, only provide upper-limits to the expected emission due to the decay of 56 Ni. These 
upper-limits on the gamma-ray flux are of 7.1 x 10 -5 ph/s/cm 2 for the 158 keV line and of 2.3 x 10~ 4 ph/s/cm 2 for 
the 812 keV line. These bounds allow to reject at the 2a level explosions involving a massive white dwarf, ~ 1 MQ in 
the sub-Chandrasekhar scenario and specifically all models that would have substantial amounts of radioactive 56 Ni in 
the outer layers of the exploding star responsible of the SN2011fe event. The optical light curve obtained with the OMC 
camera also suggests that SN2011fe was the outcome of the explosion, possibly a delayed detonation although other 
models are possible, of a CO white dwarf that synthesized ~ 0.55 Mq of 56 Ni. For this specific model, INTEGRAL 
would have only been able to detect this early 7-ray emission if the supernova had occurred at a distance & 2 Mpc. 
Conclusions. The detection of the early 7-ray emission of 56 Ni is difficult and it can only be achieved with INTEGRAL 
if the distance of the event is close enough. The exact distance depends on the specific SNIa subtype. The broadness 
and rapid rise of the lines are probably at the origin of such difficulty. 
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1. Introduction 

From the photometric point of view, Type la supernovae 
(SNIa) are characterized by a sudden rise and decay of 
their luminosity, followed by a slowly- fading tail. From 
the spectroscopic point of view, they are characterized by 
the lack of H-lines and the presence of Si II-lines in their 
spectra during the maximum light and by the presence of 
Fe emission features during the nebular phase. A noticeable 
property is the spectrophotometric homogeneity of the dif- 
ferent outbursts. Furthermore, in contrast with the other 
supernova types, they appear in all kind of galaxies. These 
properties point out an exploding object that is compact, 
free of hydrogen, that can be activated on short and long 
time scales, and is able to synthesize a minimum of 0.3 M Q 
of radioactive 56 Ni to power the light curve. These con- 
straints immediately led to the proposal that SNIa were 
the outcome of the thermonuclear explosion of a mass ac- 
cretin g C/O white dwarf (W D) near the Chandrasekhar's 
limit (|Hovle fe Fowlerl ll960) in a close binary system. 

Despite this homogeneity, when SNIa are observed in 
detail some differences appear. Now it is known that there 
is a group of SNIa with light curves showing very bright 
and broad peaks, the SN1991T class, that represents 9% of 
all the events. There is another group with a much dimmer 
and narrower peak and that lacks of the characteristic sec- 
ondary peak in the infrared, the SN1991bg class, that rep- 
resents 15% of all the events. To these categories it has been 
recently added a new one that contains very peculiar su- 
pernovae, the SN2002cx class, representing ~ 5% of the to- 
tal. These supernovae are characterized by high ionization 
spectral features in the pre-maximum, like the SN1991T 
class, a very low luminosity, and the lack of a secondary 
maximum in the infrared, like the SN1991bg class. The re- 
maining ones, which amount to ~ 70%, pr esent normal be - 
haviors and are known as Branch-normal (|Li et al.ll2011c|) . 
However, even the normal ones are not completely homo- 
geneous and show different luminosities at maximum an d 
light curves with different decline rates f|Li et al.l 1201 lblh 
This variety has recently increased with the discovery of 
SN2001ay, wh ich is charac t erized by a fast rise and a very 
slow decline (Ba ron et al.l |2012[ ). This diversity strongly 
suggests that different scenarios and burning mechanisms 
could be operating in the explosion. 

From the point of view of the explosion mecha- 
nism as seen in one dimen sional models, it is possi- 
ble to distinguish four cases (jHoeflich fe Khokhlovl fl99l 
iHillebrandt fe Niemeverll2000D : the pure detonation model 
(DET), the pure deflagration model (DEF), the delayed 
detonation model (DDT), and the pulsating detonation 
model (PDD). The equivalent models in three dimensions 
also exist, but with a larger variety of possibilities. An 
additional class are the so called Sub-Chandrasekhar's 
(SCh) models in which a detonation triggered by the 
ignition of He near the base of a freshly accreted he- 
lium layer completely burns the white dwarf. At present, 
there is no basic argument to reject any of the mod- 
els, except the DET ones that are incompatible with the 
properties of the spectrum of SNIa at maximum light. 



Present observations also pose severe constraints to the to- 
tal amount of 56 Ni that can be produced by the He-layer 
in SCh models (iHoeflich fe Khokhlovl |1996t iNugent et all 
119971: IWooslev fc Kasenll2011l) . 

According to the nature of the companion, either 
non-degenerate or degenerate, progenito rs can be clas- 
sified as single degenerate systems -SD (IWhelan fe Ibenl 
1973) or double degen erate systems -DD (|WebbinkHl984t 
llben fe Tutukov| [l985). The distinction among them is im- 
portant in order to interpret the observations since, de- 
pending on the case, the white dwarf can ignite below, near 
or above the Chandrasekhar's mass and consequently the 
total mass ejected and the mass of 56 Ni synthesized can 
be different. At present, it is not known if both scenarios 
can coexist or just one is enough to account for the su- 
pernova variety. Observations of the stellar content in the 
interior of known SNIa remnant s point towards one possible 
SD candidate (Tycho SNR see Ruiz-Lapuente et ahl f2004: 
iGonzalez Hernandez et all 120091 : iKerzendorf et al.l 12009) 
and two almost certain DD candidates (SNR0509-67.5 and 
SNR 0519-69.0 ISchaefer fe Pagnottal I20TI lEdwards et"ai1 
I2012h . 

The detection of 7-rays from supernovae can provide 
important insight on the nature of the progenitor and es- 
pecially on the explosion mechanism, since the amount and 
distribution of the radioactive material produced in the 
explosion strongly depend on how th e ignition starts and 
how the nuclear flame propagates (see lGomez-Gomar et al.l 
119981 : llsern et al.l2008l for a detailed discussion of how these 
differences are reflected in the spectra). The advantages of 
using 7-rays for diagnostic purposes relies on their penetra- 
tive capabilities, on their ability to distinguish among differ- 
ent isotopes and on the relative simplicity of their transport 
modelling as compared with other regions of the electro- 
magnetic spectrum. Unfortunately, such observations have 
not bee achieved so far because of the poor sensitivity of 
the instruments. For this reason up to now it h as only been 
possi ble to place uppe r limits to the SN19 91T ([Lichti et al.l 
I1994D and SN1998bu (jGeorgh et al.ll2002D events. 

Several authors have exami n ed the 7-ray emission 
from S NIa dGehrels et al.l 119871: lAmbwani fe Sutherlanc 
| 1988|: iBurrows fe Thej 09m iRuiz-Lapuente etal 



Kumagai & Nomotol 11997 



* Based on observations with INTEGRAL, an ESA project 
with instruments and science data centre funded by ESA 
member states (specially the PI countries: Denmark, France, 
Germany, Italy, Switzerland, and Spain), the Czech Republic, 
and Poland and with the participation of Russia and USA. 



1 19931: iHoeflich et all 11994 

Tim mes fe Wooslevl Il997t iGomez-Gomar et al.l 11998 
ISim fe Mazzalill2008l) . To explore the above model variants 
we have used as a g uide the properties o btain ed with 
the code described in IGomez-Gomar et al.l (11998 1 ), which 
i s bas ed on the methods described by Pozdnvakov et all 
(119831 ) and lAmbwani fe Sutherland] (|l988l h In order to 
test the consistency of this model, the results of this code 
were successfull y cross-checked w ith those obtained by 
other authors ( Milne et al.l l2004|). This code w as later 
generalized to three dimensions ~ pirschmanr]|2009D . 

Before and around the epoch of maximum of the optical 
light curve, the 7-ray emission can be characterized (figure 
2 of Gomez-Gomar et al. loc.cit.) as follows: i) A spectrum 
dominated by the 56 Ni 158 and 812 keV lines, ii) Because of 
the rapid expansion, the lines are blueshifted but their en- 
ergy peak quickly evolves back to the red as matter becomes 
more and more transparent. The emergent lines are broad, 
typically from 3% to 5%. Because of the Doppler effect the 
812 keV line blends with the quickly growing 56 Co 847 keV 
line, forming a broad feature, iii) The intensity of the 56 Ni 
lines rises very quickly with time, being very weak at the 
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beginning, even in the case of Sub-Chandrasekhar models. 
This fact, together with the relatively short lifetime of 56 Ni, 
makes the observational window rather short. 

SN 2011fe (RA = 14:03:05.81, Dec = +54:16:25.4; 
J2000) was discover ed in M101 on August 24th, 2011 
(| Nugent et al.l l2011aT) . The absence of hydrogen and he- 
lium, coupled with the presence of silicon in the spectrum 
clearly indicates that it belongs to the SNIa class. Since it 
was not visible on August 23th, this supern ova must have 
been d etected ~ 1 day after the explosion (|Nugent et al.l 
l2011al). Furthermore, as M101 is at a distanc e of 6.4 Mpc 
(|Stetson et aII[T99l IShappee fc Stanekll2011ft . SN2011fe is 
the brightest SNIa detected in the last 25 years. This dis- 
tance is slightly less than the maximum distance at which 
current gamma-ray instruments should be able to detect an 
intrinsically luminous SNIa. The closeness of SN2011fe has 
made it possible to obtain the tightest constraints on the 
supernova and its progenitor system to date in a variety of 
observational windows. Red giant and helium stars compan- 
ions, symbiotic systems, systems at the origin of optically 
thick wind s or containing recurrent novae are excluded for 
SN20 11fe (iLi et alj|2011at iBloom et alj|20ll iB rown et al.l 
l201lt IChomiuk et alJ 1 12012} ). leaving only either DD or a 
few cases of SD as possible progenitor systems of this su- 
pernova. 

In this study we analyze the data obtained by 
INTEGRAL during the optical pre-maximum observations 
spanning from 4258.8733 IJD (August 29th, 20:59 UT) to 
4272.1197 IJD (September 12th, 2:52 UT) with a total ob- 
servation time of 975419 s. This schedule was essentially 
determined by the constraints imposed by the Sun, accord- 
ing to the TVP tool of INTEGRAL, that prevented the ob- 
servation just beyond the optical maximum where the 56 Ni 
lines are expected to peak. Despite this limitation, these 
early observations were triggered to constrain any predicted 
early gamma-ray emission as may be expected from some 
variants of sub-Chandrasekhar models. In the next section, 
we present the data obtained with the instruments onboard 
INTEGRAL. Then, we discuss the limits they can put on 
present models of SNIa explosion, and conclude. 



2. INTEGRAL data 



INTEGRAL ()Winkler et al.1 12003D is able to observe in 
gamma-rays, X-rays and visible light. It was launched in 
October 17th 2002 and was injected into a highly eccentric 
orbit with a period of about 3 days in such a way that it 
spents most of its time well outside the radiation belts of the 
Earth. The spacecraft contains two main instruments, SPI, 
a germanium spectrometer for the energy range of 18 keV 
to 8 MeV with a spect ral resolution of 2.2 keV at 1.33 MeV 
(jVedrenne et al.ll2003ft . and IBIS, an imager able to provide 
an im aging resolution of 12 arcmin FWHM (jUbertini et al.1 
12003ft . which has a CdTe detector, ISGRI, able to provide 
spectral information of the co ntinuum and broad lines in 
the range of 15 keV to 1 MeV ()Lebrun et alJl2003t ). Other 
instruments onboard are an X-r ay monitor, JEM -X, that 
works in the range of 3 to 35 keV (jLund et al.ll2003ft . and an 
optical camera, OMC, able to operate in the visible band 
of th e spectrum up to a magnitude of 18 (|Mas-Hesse et al.1 
12003ft . 
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Fig. 1. Light curves in the V-band obtained from delayed 
detonation models that produce, up to down, 0.63, 0.45, 
0.27, 0.09 M of 56 Ni and satisfy the brightness-decline 
relationship. 



2.1. OMC data 

The height of the maximum of the optical light curve de- 
pends on the total amount of 56 Ni synthesized during the 
explosion, its distribution in the debris, the to tal kinetic en- 
ergy of matter and the opacity (|Arnettlll996[ ). The bright- 
ness decline relation ( Amis) relates the absolute brightness 
at maximum light and the rate of the post-maximum de- 
cline over 15 days. From theory, this relationship is well un- 
derstoo d: light curves are powere d by the radioactive decay 
of 56 Ni (jCokate fc McKeel 1 1969ft . More 56 Ni increases the 
luminosity and causes the envelopes to be hotter. Higher 
temperature means higher opacity and, thus, longer dif- 
fusion tmiB_s£ales_and_slower_o^ 

light (iHoefhch et al.ll996tlNugent et alJll997tlUmeda et al.l 
1999: iKasen et afll2009ft . The Ami5-relation holds up for 
virtually all explosion scenarios as long as there is an excess 
amount of stored e nergy to be released (|Hoeflich et al.ll995 
iBaron et al.1 [20 1 2 i. The tightness of the rela tion observed 
for Branch- normal SNe la is about 0.3 m (Ha muv et al.l 
119961 : iPerlmutter et al.1 11999ft and it is consistent with ex- 
plosions of models of similar mass. Since delayed detonation 
models (see Fig. [lj provide a reasonable spectral evolution, 
their maximum brightness is consistent with the Hubble 
constant and the Am 15 -relationship is consistent with the 
observations with a factor between 1 and 1.3, we will use 
these models to estimate the total amount of 56 Ni freshly 
synthesized. 
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Fig. 2. Light curve of SN2011fe in the V-band obtained 
with the OMC / 'INTEGRAL assuming a distance of 6.4 Mpc 
(dots) for the orbits 1086-1088 and 1097-1101. The contin- 
uous line represents a reasonable fit obtained with a DDT 
model (see text). 



The properties and t he photometric characte rization of 
the OMC can be found in lMas-Hesse et al.l (|2003h . The data 
obtained during the orbits 1084-1088 plus the data corre- 
sponding to orbits 1097-1111 were analysed with the Offline 
Scientific Analysis Software (OSA, versi on 9) provided by 
the I SDC Data centre for Astrophysics (|Courvoisier et al.l 
120031) . The fluxes and magnitudes were derived from a pho- 
tometric aperture of 3 x 3 pixels (1 pixel = 17.504 arcsec), 
slightly circularized, i.e. removing 1/4 pixel from each cor- 
ner (standard output from OSA). The default centroiding 
algorithm was used, i.e. the photometric aperture was cen- 
tred at the source coordinates. We checked that the photo- 
metric aperture of 3 x 3 pixels does not include any signifi- 
cant contribution by other objects. Because the source was 
bright enough, combination of several shots were not re- 
quired to increase the signal-to-noise ratio. To only include 
high-quality data, some selection criteria were applied to in- 
dividual photometric points removing those measurements 
with a problems flag. Shots were checked against satura- 
tion, rejecting those with long exposures (200 seconds) for 
V < 10.5. To avoid noisy measurements the shortest expo- 
sures (10 seconds) were not used. 

The extinction along the line of sight has two compo- 
nents, one due to the Milky Way and other to M101. Since 
M101 has a galactic latitude of ~ 60° the Milky Way con - 
tribution is expected to be small. IShappee fe Staneld ( 201 it ) 
estimate E(B — V) = 0.009 which means an extinction of 
Ay ~ 0.03 mag. The M101 contribution is also expected 



to be small since it is a face on galaxy and SN2011fe is 
placed in a regio n with a small concentra tion of interstel- 
lar dust and gas ([Suzuki et al.ll2007l I2009T) . It is possible to 
obtain a rough upper limit of the extinction using the obser- 
vations (Sha ppee fe Stanekll20lil ) of two regions containing 
Ccpheids that overlap just at the position of SN2011fe. The 
average of the two extinctions, E(B — V) = 0.2 mag, would 
imply an absorption of Ay = 0.62 mag. However, since the 
region containing the supernova has a small concentration 
of interstellar matter, the comparison of the light curve at 
different bands and the absence of strong sodium lines in 
the spectrum suggest that extinction is small. Thus the 
adopted absorption affecting the supernova is taken to be 
Ay ~ 0.03. 

The V-band light curve reached the maximum, V — 
9.99, at the day IJD = 4271.44, in agreement with 
iRichmond fe Smith! (l2012f). Taking a canoni c al di stance 
of 6.4 Mpc (see however Tammann fc Reindll (120111 )) and 
assuming no extinction, the absolute magnitude should 
be My = —19.04 at maximum thus indicating that the 
SN2011fe was a slightly dim average SNIa. The V-band 
light curve can be well fitted with a delayed detonation 
model (see Fig. [5]) of a Chandrasekhar mass WD igniting 
at pc = 2 x 10 9 gcm -3 , making the transition deflagra- 
tion/detonation at p tr — 2.2 x 10 7 gcm -3 . This model pro- 
duces 0.51 M Q of 56 Ni, although if extinction and distance 
uncertainties are taken into account this v alue could eas- 
ily be a - 10% larger (jrloflich et alJl2002ft . This value is 
also c onsistent with the estimations found bv lRopke et al.l 
(2012) using an independent DDT model and a violent 
merger model and is roughly equivalent to the one obtained 
with model DDTe of Table d We note that this yield of 
56 Ni is consistent with the value derived by iNueent et al.l 
(|2011bD . From this theor etical model, Am 15 ( £) = 1.2 ± 



0.2, in agreement with iTammann fc Reindll (12011) but 
slightl y smaller than the value found bv lRichmond fc Smith! 
(|2012l). The observ ation of the e arly light curve, 4 ho urs 
([Bloom et al.ll2012l ). and 11 hours (|Nugent et alJl2011bD af- 
ter the explosion strongly supports the scenario based on 
the explosion of a C/O white dwarf. 

2.2. SPI data 

Only data of the low energy range (<^ 2 MeV) have been 
analysed in this study of the SPI output. The energy cali- 
bration was performed for each orbit by fitting parameters 
of a four degree polynomial function with the channel po- 
sitions of the instrumental background lines at 23.4 keV, 
198.4 keV, 309.9 keV, 584.5 keV, 882.5 keV and 1764.4 keV. 
The precision of the resulting calibration is better than 
~ 0.1 keV at 1 MeV. The calibrated single-detector and 
multiple-detectoiQevent data have been binned into sepa- 
rated spectra at 0.5 to 50 keV per bin. 

The time-averaged energy spectrum of the supernova 
has been extracted by a model fitting method. The flux of 
the source and the instrumental background are both fitted 
to the data (counts per detector per pointing at each en- 
ergy bin), assuming a point source at the SN2011fe position. 
The instrumental background is fitted per pointing assum- 
ing that the relative background rate in each detector of the 



1 We used only double-detector events as including events with 
larger detect or multiplicity doe s not improve the sensitivity for 



this analysis (jAttie et all 2003). 
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Fig. 3. Time variation of the event rate in the ACS of SPI 
during the pre-maximum observations of SN2011fe. 
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Fig. 4. Time variation of the event rate in the GeD of SPI 
during the pre-maximum observations of SN2011fe. 



Ge camera is fixed, and is obtained by summing the counts 
per detector of all the pointings of the observation. We have 
verified that the detector pattern did not change with time. 
This background modelling method is adapted to the anal- 
ysis of data that show strong instrumental background vari- 
ations with time scales ^ 3 dayfQ. Figures [3] and E] display 
the count rate in the anticoincidence system (ACS) and in 
the germanium detector (GeD) of SPI. The ACS count rate 
could be used to trace the instrume ntal background fl uctu- 
ations in the germanium detectors ()Jean et~atl[200l . The 
passage through the radiation belts and the impact of the 
solar activity are clearly seen. Solar activity was particu- 
larly influential during this period of observations because 
of the occurrence of two solar flares on ~ 4266 and 4267 
IJD followed by a Forbusch decrease around 4268.5 IJD. 
In order to check that these background variations do not 
produce artifacts, we also performed the model fitting anal- 
ysis using a filtered data set without periods with strong 
rate variations. The results obtained with filtered and non- 
filtered data are statistically equivalent. Consequently, we 
decided to use non-filtered data in the next steps of the 
analysis. 

Figure [5] displays the spectrum of the whole observation. 
It has been obtained by combining the spectra extracted 
by model fitting from single-detector and multiple-detector 
events. In this case, the model fitting has been performed 
with data rebinned in 50 keV bins. The spectrum does not 
show any significant feature. A \ 2 test shows that it is 
consistent with a Poissonian background (x 2 = 12.4 with a 
dof = 10). A similar conclusion is obtained when the spectra 
are extracted by model fitting with data rebinned in 5 keV 
or 2 keV bins. 

No significant excess was found in the spectrum, even 
at the energies of the strongest gamma-ray lines that are 
expected from the decay of 56 Ni (blueshifted 158 and 812 
keV lines) . Table [T] presents the upper- limit fluxes derived 
from th e analysis, at the e nergies of interest for several band 
widths (jlsern et al.ll2011l) . 



The analyses performed with background models that use 
background tracers show strong systematics due to these large 
instrumental background variations. 
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Fig. 5. Spectrum of the whole early observation of 
SN2011fe extracted from SPI data with 50 keV bins. 

Table 1. Upper- limit of the flux in selected spectral regions 
for SPI (2a), JEM-X (2a), and IBIS/ISGRI (3a) for the 
entire pre-maximum observation period. 

Energy band Upper-limit flux Instrument 
(keV) (photons s _1 cm~ 2 ) 
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2.3. IBIS/ISGRI data 

IBIS uses two detection layers to cover the same energy 
range as SPI. The low energy camera, ISGRI, uses 16,384 
thin CdTe detectors operated at ambient temperature. The 
imaging is much better (12' FWHM) than for SPI but the 
spectral resolution is more limited and the efficiency begins 
to drop above 100 keV. As a result, the IBIS/ISGRI sensi- 
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tivity to the 812 keV line from 56 Ni is far worse than the 
SPI one at the same energy even in the case of a broad line. 
Conversely, IBIS/ISGRI is much better at work than SPI 
to reveal a low energy continuum or broad lines such as 
the 158 keV 56 Ni or 122 keV 57 Co ones. To summarize the 
picture, one could say that SPI is better at hand to detect 
narrow lines in a spectrum while IBIS/ISGRI is better at 
detecting point sources in broad energy-range sky images. 
Two data processing methods have been followed in par- 
allel. One uses the standard OSA-9 version and the other 
takes advantage of the developments for the forthcoming 
OSA-10. The latter approach includes two new corrections: 
for the spectral drift along the mission and for the flat field. 
In either case no significant (greater than 3 a) signal was 
found at the position of SN2011fe. Table [T] gives the up- 
per limits obtained for the 60-170 keV, 90-172 keV and the 
150-172 keV bands. Similarly, maps using the last 6 days 
of the observing period for the same energy ranges do not 
show any significant emission at the position of SN2011fe. 

2.4. JEM-Xdata 

Both JEM-X units were simultaneously operating at the 
time of these INTEGRAL observations and they can be 
used to constrainthe continuum emission using a source 
search in broad band images as in the IBIS/ISGRI case. 
Due to the smaller field of view of the JEM-X mon itors 
compared to those of SPI and IBIS dLund et al.ll2003l ). we 
selected and analyzed only those pointings where SN2011fe 
was within 5° of the pointing direction. The data obtained 
during the orbits 1084 - 1088 (August 29th to September 
12th, 2011) were analysed with OSA-9 following the stan- 
dard procedure. Images from single pointings were com- 
bined into one mosaicked image for each X-ray monitor. 
The two composite images from JEM-X1 and JEM-X2 
were then merged to obtain the final image, providing an 
on-source effective exposure time of 450 ks. The imag- 
ing analysis was performed in the 3-25 keV band, and 
in the 3-10 keV and 10-25 keV sub-bands. SN2011fe is 
not detected in any of the JEM-X images. Assuming a 
Crab-like spectrum, we estimate 2 o upper limits on the 3- 
10 keV and 10-25 keV fluxes of 5 x lO^phcm- 2 s" 1 and 
4x 10 _4 phcm -2 s , respectively (see Table[TJ). The source 
is also not detected when analyzing separately the JEM-X 
data from orbits 1086 - 1088 (Sept 4th to 12th, 2011), cor- 
responding to the observations closer to the epoch of the 
expected 56 Ni line maximum. 

3. Discussion 

The 7-ray spectrum of SNIa mainly depends on the total 
amount and distribution of 56 Ni synthesized during the ex- 
plosion, as well as on the chemical structure and ve locity 
distribution of the debris (|G6mez- Gomar et all [19981. For 
instance, when 56 Ni is present in the outer layers of some 
of the sub-Chandrasekhar models, the corresponding lines 
should appear very early in the spectrum. However, for a 
given flux, there are at least two more factors that deter- 
mine if the 7-signal is detectable: the change of the signal 
with time and the width of the lines. 

In order to check the influence of the mass and dis- 
tribution of 56 Ni, several models obtained under differ- 
ent hypothesis about the burning regime or the explo- 
sion mechanism have been considered. Most are one- 



Table 2. Kinetic energy (K) and mass of 56 Ni produced by 
different models of explosion. 



Model 


K (foe) 


Mm (Mq) 


DETO 


1.44 


1.16 


DD202c 


1.30 


0.78 


DDTc 


1.16 


0.74 


SC3F 


1.17 


0.69 


W7 


1.24 


0.59 


SCOP3D 


1.17 


0.56 


DDTe 


1.09 


0.51 


SC1F 


1.04 


0.43 


HED6 


0.72 


0.26 



dimensional spherically symmetric models, which should 
be sufficient for SN2011fe in view of the small amount of 
global asy mmetry suggested by spectropolarimetric mea- 
surements (|Smith et al.ll201ll ). We have not i ncluded mod- 
els of D P explosions for two reasons. First, iNueent et al.l 
(201 lbj) reproduced satisfactorily the rising part of the light 
curve of SN2011fe with a simple analytic model involving 
a Chandrasekhar mass WD, and found only small amounts 
of unburned carbon in the early spectra of this supernova. 
Second, in spite of recent advances on simula tions of almost 
norm al SNIa from mergers of massive WDs (jPakmor et al.l 
l2012f ). theoretical models of DD explosions are not as ma- 
ture at present as those involving either Chandrasekhar 
or sub-Chandrasekhar mass WDs accreting from a non- 
degenerate companion. Hence, the number of free parame- 
ters involved in DD explosions is too large to allow for effi- 
cient constraints derived from upper limits on the gamma- 
ray emission of SN2011fe. 

Table [2] displays the main characteristic s of the models 
used in the present study. The DETO model (jBadenes et al.l 
2003) corresponds to a pure detonation of a WD near the 
Chandrasekhar' s limit. It is al so representative o f the most 
massive models compu ted bv iFink et all (l201fl . The W7 
i s the classical mode l of lNomoto et all (ll984~The DD202c 
(jHoeflich et al.lll998D and the DDTc,e (jBadenes et al.12 0051 
models are delayed detonation models that produce dif- 
ferent amounts of 56 Ni and have different expansion en- 
ergies. The HED6 model corresponds to the explosion of 
a 0.6 Mfn C/O WD that has acc reted 0.17 M Q of he- 
lium of iHoeflich fc Khokhlovl (fl996lh and SC1F and SC3F 
(E. Bravo, unpublished) are also sub-C handrasekh a r mod - 
els equivalent to models 1 and 3 of IFink et ahl ((2010). 
Finally, the SCOP3D model is a three dimensional sub- 
Chandrasekhar model t hat corresponds to model A of 
iGarcia-Senz et al.l (|l999h . 

Figure [B] displays the SPI 2<r upper-limit spectrum ob- 
tained with 5 keV bins as well as the gamma ray spectra 
predicted by several models for this early observation, as- 
suming that the distance of SN2011fe is 6.4 Mpc. All the 
models are well below the upper-limit. The predicted in- 
tensity of the 56 Ni 158 keV, 270 keV, 480 keV, 750 keV, 
and 812 keV lines is maximum for the detonation model 
(DETO) and the Sub-Chandrasekhar's models SCOP3D, 
SC3F, and SC1F, as expected. Since the spectral shape 
(width and centroid) and intensity change from one model 
to another, the energy band used to extract the fluxes of 
every line or complex of lines was chosen to provide the 
optimum significance for each model. Results are displayed 
in Table El 
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Table 3. Measured and predicted flux in spectral regions where 56 Ni lines are expected to be as a function of models. 
The fifth column is the probability to reject the model (see text). In the case of DDTc, W7 and DDTe models, the fluxes 
were obtained for an optimal observation periods of 7, 5 and 5 days, respectively (see Tabl^31). 



Model 


Energy band 


Measured flux 


Predicted flux 


Probability 




(keV) 


(10~ 4 photons s _1 cm -2 ) 


(10 -4 photons s~ x cm -2 ) 


% 


DETO 


150 - 


175 


-0.20 ± 0.78 


0.81 


92.4 




265 - 


295 


-0.55 ± 0.85 


0.48 






482 - 


508 


0.13 ± 0.88 


0.36 






730 - 


880 


-1.01 ± 2.02 


2.01 




SCOP3D 


145 - 


170 


-0.50 ± 0.76 


0.70 


84.7 




267 - 


297 


-0.24 ± 0.86 


0.36 






482 - 


508 


0.13 ± 0.88 


0.27 






720 - 


880 


-0.75 ± 2.08 


1.54 




SC3F 


150 - 


170 


-0.03 ± 0.67 


0.37 


62.4 




253 - 


283 


-0.99 ± 0.85 


0.25 






452 - 


508 


0.38 ± 1.42 


0.29 






720 - 


880 


-0.75 ± 2.08 


0.76 




SC1F 


150 - 


170 


-0.03 ± 0.67 


0.32 


51.8 




254 - 


288 


-0.81 ± 0.90 


0.22 






458 - 


509 


0.01 ± 1.37 


0.23 






725 - 


900 


-1.02 ± 2.18 


0.63 




HED6 


159 - 


169 


-0.13 ± 0.49 


0.08 


21.3 




270 - 


284 


-0.98 ± 0.59 


0.05 






475 - 


505 


-0.27 ± 1.04 


0.06 






746 - 


872 


-0.13 ± 1.87 


0.27 




DDTc 


158 - 


165 


-0.10 ± 0.57 


0.03 


6.2 




270 - 


282 


-0.58 ± 0.76 


0.02 






480 - 


498 


-1.22 ± 1.00 


0.02 






740 - 


880 


1.28 ± 2.08 


0.08 




W7 


158 - 


165 


-0.34 ± 0.69 


0.02 


3.4 




270 - 


282 


-0.88 ± 0.92 


0.02 






478 - 


498 


-1.82 ± 1.37 


0.02 






740 - 


880 


-0.31 ± 1.38 


0.04 




DDTe 


158 - 


165 


-0.34 ± 0.69 


0.01 


2.0 




270 - 


280 


-0.49 ± 0.84 


0.008 






480 - 


498 


-1.00 ± 1.26 


0.009 






740 - 


880 


-0.152 ± 3.25 


0.05 





Despite the non-detection and the bounds imposed by 
the optical light curve, the compatibility of the different 
models with the zero flux hypothesis has been tested. Since 
the measured fluxes are compatible with zero, the x 2 value 
for each model was computed as the quadratic sum of the 
expected significance of the predicted fluxefl Then, the 
probability of rejection of each particular model was de- 
rived by computing the probability that the % 2 is smaller 
than the measured value within \ 2 statistics. The rejection 
probability of the selected models is given in the fifth col- 
umn of Table [3] None of them can be firmly rejected by SPI 
data. Although the probability to reject the DETO model 
is 92.4% this value corresponds to a significance of ~1.4cr. 
To reject it at a 99% of confidence level, the intensity of 
the lines should be larger by a factor 1.45 times or, equiv- 
alently, the distance to SN2011fe should be less than 5.3 
Mpc. Since the DETO model is the one that produces the 
brightest lines, SPI can only provide interesting results dur- 
ing this epoch if the distance to the supernova is smaller 
than this value. 

The influence of the temporal behavior of the line inten- 
sity and line width on the detectability of the 7 emission 
can be easily seen by estimating the significance of the ob- 



The expected significance is the predicted flux divided by 
the uncertainty of the flux measurement in the corresponding 
energy band. 



servation. In the lim it of weak signals this significance is 
given by (|Jeanlll996lh 



A eff J <p(t)dt 



VbVAEAt 

where At is the observation time, A e ff the effective area at 
the corresponding energies, ip is the flux (cm _2 s _1 ) in the 
energy band AE, V is the volume of the detector and b is 
the specific noise rate (cm _3 s~ 1 keV _1 ), where it has been 
assumed that it is weakly dependent of the energy and time 
in the interval of interest. 

If the flux grows like <p (t) = ifioe at , the significance 
reached by integrating the time interval (U — At,ti) is 



-aAt 



A e f f (p(tj) l-e 
VabVAE VaAl 



(2) 



For aAt << 1, it behaves as n oc v At and has a maxi- 
mum at aAt = 1.26. Furthermore, the dependence on AE, 
clearly shows the convenience of taking a value that maxi- 
mizes the signal/noise ratio. Unfortunately, since the value 
of a is not known a priori the optimal observing time is not 
known in advance. Since the time dependence of the model 
fluxes does not follow strictly an exponential growth, we 
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Fig. 6. Comparison of the mean flux obtained from the the- 
oretical models during the early observation period with the 
2cr upper-limit flux in 5 keV bins. Th eoretical spectra were 
obtain ed with the code described in iGomez-Gomar et al.1 



have analysed the optimum observation periods for each 
model. Table S] displays the results. In general, models and 
energy ranges showing weak variation of the flux are best 
detected when the observation period is large whereas mod- 
els with a strong variation of their flux during the first two 
weeks are best detect ed using data at those time s. 

It is well known (Go mez-Gomar et al.1 Il998l ) that the 
width of the line has a strong influence on their detectabil- 
ity. In the case of the present observations it is possible to 
estimate the significance of a narrow line, AE ^ 2 keV, 
by comparing the measured la flux uncertainty (see Table 
[3]), ip i, with the flux ip and the width AE predicted by the 
model 



AE 



tpi V 2 keV 



(3) 



In the case of the 158 keV line (Tabled, the flux predicted 
by the DETO model is <p = 8.1 x 1CT 5 ph s^cm" 2 , ipx = 
7.8 x 1CT 5 ph s _1 cm~ 2 , and AE = 25 keV, where AE has 
been chosen to maximize the signal/noise ratio. With these 
data, n ~ 3.8 and in this case, with this hypothesis, the 158 
keV line would have been detected by SPI if it had been 
narrow. On the contrary, in the model SC1F, ip = 3.2 x 10~ 5 
ph s-^m" 2 , tpt = 6.7 x 1CT 5 ph s^cm -2 , AE = 20 keV 



Table 4. Optimal observation period for several models 
and energy ranges with SPI. 



Model 


Energy band 
(keV) 


Optimal period 
(days) 


DETO 


150 - 175 


13.5 




730 - 880 


13.5 


SCOP3D 


145 - 170 


13.5 




750 - 880 


13.5 


SC1F 


150 - 170 


13.4 




725 - 900 


13.4 


HED6 


155 - 175 


11.7 




730 - 880 


12.4 


DDTc 


70 - 165 


8.0 




740 - 880 


6.2 


W7 


70 - 900 


5.7 




820 - 840 


4.7 


DDTe 


158 - 165 


4.9 




740 - 880 


4.9 



7 10" 

% io- 



m 10- 



io- 



10" 





III 

fey 




1 ■ 







0.5 



1.5 



Energy (MeV) 




0.75 



0.8 0.85 
Energy (MeV) 



0.9 



Fig. 7. Early gamma ray spectra predicted by the delayed 
detonation model that best fits the OMC (see Fig. at 
days 8, 11, 14, 17, 20 (green, blue, cyan, magenta, black) 
at a distance of 6.4 Mpc (upperer figure). The lower pannel 
displays the evolution of the profile of the 56 Ni 812 kev- 
56 Co 847 keV features. 



and the result is n ~ 1.5. Here, the line would only be 
marginally detectable even in the narrow case. 

The luminosity at the maximum of the optical light 
curve is proportional to the total mass of 56 Ni, while the 
shape of the optical peak is determined by the amount o f 
56 Ni and the kinetic energy of the remnant (|Arnettlll997ft . 
As it has already been mentioned, the favoured model by 
the optical light curve is the DDTe one, and it seems worth- 
while to interpret the SPI observations in terms of this 
model. 

Figure [7] (upper panel) displays the early gamma ray 
spectra predicted by DDTe for several instants after the 
explosion. The way the different lines grow and the width 
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of them is clearly seen. The main problem with the 812 
keV line is that, due to the Doppler effect and energy 
degradation of photons by Compton scattering, it blends 
with the 847 keV line of 56 Co (see Fig. [71 lower panel) 
and makes more difficult the interpretation of the signal. 
During this observation, and for this model, the energy in- 
terval that maximizes the signal/noise ratio is ~ 25 keV 
centered at 826 keV, and the corresponding flux at day 18 
is <p = 3.4 x 10~ 6 ph s^'crn" 2 . On the contrary, the 158 
keV line has a more regular behavior and it is better suited 
for diagnostic purposes. At day 20 after the explosion, the 
flux in this line is (p = 4.7 x 1CT 6 ph s"W 2 for AE = 20 
keV. In both cases, the signal is too weak to be detected 
with SPI and the source should be at a distance of ^ 2 Mpc 
to be detectable. 

The 200-540 keV band contains almost all the photons 
produced by the annihilation of positrons produced during 
the 56 Ni chain decay. Since the emission of 511 keV only 
represents the 25% of the total number of the annihilation 
photons, this band offers, in principle, better possibilities of 
detection than the 511 keV line itself. In the DDTe model, 
the maximum of the band occurs at day 50 and amounts 
1.05 x 10~ 4 ph s _1 cm~ 2 , while the maximum of the 511 
keV line occurs at day 74 with a flux of 1.2 x 10 -5 ph 
s _1 cm -2 , within a band of 30 keV. The expected emission 
at day 18 after the explosion, the last day of this observa- 
tion window, are 4.9 x 10 -5 ph s _1 cm~ 2 and 2.4 x 10~ 6 
ph s _1 cm -2 respectively. The 3cr sensitivity in this band is 
estimated to be 8.7 x 10 -4 ph s _1 cm -2 for an integration 
time of 10 6 s, well above the value of the expected emission 
for this model, for which reason this band is not detectable 
by SPI/ INTEGRAL except if the explosion occurs at a dis- 
tance d Si 1 — 2 Mpc (assuming an explosion similar to that 
predicted by the DDTe model). 

As in the case of SPI, it is possible to discuss the results 
from ISGRI in more depth using detailed predictions by 
different theoretical models. In this case we use the model 
to produce a daily spectrum corresponding to the period 
of observation and combine these spectra with the ISGRI 
sensitive area (ARF) and spectral response (RMF) to pre- 
dict the observable ISGRI spectrum as a function of time. 
Using the latest ISGRI background estimate (F. Mattana, 
private communication), the expected relative value of sig- 
nal to noise ratio can be computed as a function of energy 
and time. This was done for the DDTe, W7, DETO and 
SC3F models. The first two of these models predict fluxes 
that are out of reach of the ISGRI sensitivity by an order 
of magnitude even if only the last days of the observing pe- 
riod are used, where the expected signal is higher. Figure 
\E\ illustrates an attempt in the case of the DDTe model. 
The 60 - 172 keV upper limit on the ISGRI rate during the 
last six days of the first observation period is 1.7 x 10~ 3 
s _1 kcV _1 , which represents 2.2 x 10~ 4 cm~ 2 s _1 . 

Since the DETO and SC3F models produce important 
amounts of 56 Ni near the surface, they display an impor- 
tant emission feature in the ~ 150 — 170 keV band. Figure 
[5] displays the spectrum that could be expected from the 
DETO model. This spectrum is an average over all the time 
that ISGRI has been pointed to the source and the opti- 
mal band corresponds to 150 - 172 keV. In this case the 
3 a limit corresponds to 7.1 x 10~ 5 cm _2 s _1 . If a source 
like DETO would have been present, the significance of 
the signal would have been 2.8 er, but there was nothing 
detected by ISGRI in the field. Model SC3F exhibits a sim- 
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Fig. 8. ISGRI expected rate-spectrum from a DDTe model 
averaged over the last 6 days of the first observation period. 
The 3rj ISGRI upper limit for the 60-172 keV band is dis- 
played as an horizontal bar. 
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Fig. 9. ISGRI expected rate-spectrum from a DETO 
model averaged over the entire observation period. The 3 a 
ISGRI upper limit for the 150-172 keV band is displayed 
as an horizontal bar. 

ilar behavior (figure [lU] ) and, for the same specifications, 
the significance in the band is 1.4 a. However, in this case, 
the continous emission at low energies is important and the 
optimal observation band is 90 - 172 keV, the 3cr limit is 
1.1 x 10~ 4 cm _2 s^ 1 and the significance of the signal would 
have been 2.1 er. Notice that if these models, that synthe- 
size 1.16 and 0.69 M of 56 Ni respectively, are scaled to the 
~ 0.5 M Q demanded by the observations of SN2011fe in the 
optical, their expected gamma-ray flux would be marginally 
above the ISGRI sensitivity limit. 

4. Conclusions 

SN2011fe has been observed with the four instruments 
SPI, ISGRI/IBIS, JEM-X and OMC operating on board 
of INTEGRAL just before the maximum of the optical 
light curve for a period of time of 1000 ks and the data 
were compared with the predictions of several theoretical 
models. SPI data in the bands containing the 158 keV and 
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Fig. 10. ISGRI expected rate-spectrum from a SC3F 
model averaged over the entire observation period. The 3 a 
ISGRI upper limit for the 90-172 keV band is displayed as 
an horizontal bar. 

the 812 keV emission of 56 Ni would have allowed to reject 
at 98% confidence level (2er) models that produce large 
amounts of 56 Ni near the surface, if the supernova were 
closer. For instance, models with a mass fraction of 56 Ni 
in the outer 0.15 M Q of 0.05 and 0.5 would have been re- 
jected if the supernova distance was closer than 1.4 and 3.7 
Mpc respectively, assuming non-detection. Furthermore, 
ISGRI/IBIS has proved to be very efficient to explore the 
low energy region (<~ 100 keV) and has confirmed that there 
were not significant amounts of radioactive elements in the 
outer layers. This picture is consistent with the light curve 
obtained with the OMC. The observations in the optical 
suggest that the total amount of 56 Ni produced in the event 
is <~ 0.5 M originated by a mild delayed detonation explo- 
sion. 
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